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Stationary thick brane solutions supported by two spinor fields are considered. Two spinor fields 
are used here to exclude the off-diagonal components of the energy-momentum tensor of the spinor 
fields. The trapping of a test scalar field on the brane is also considered. 



I. INTRODUCTION 

At present, almost all physical theories pretending to the role of a theory unifying all known fundamental interactions 
are formulated in spacetime having more than four dimensions. Extra dimensions may have either microscopic size 
of order of the fundamental Planck length or they may be large, and even infinite. In the first case, the smallness of 
extra dimensions can be achieved by introducing special mechanisms providing compactification of extra dimensions 
when their characteristic size becomes much less than that of four-dimensional spacetime. 

In models with large extra dimensions, it is assumed that all matter fields are confined on some hypersurface (called 
a brane) which, in turn, is embedded in some multidimensional space (called a bulk) [J. In such approach, two classes 
of branes are introduced: infinitely thin branes and thick brancs. A brane is called infinitely thin if it has a delta-like 
distribution of matter along the extra dimension. However, from a realistic point of view, it is most reasonable to 
expect that a brane has a finite thickness. There are many works devoted to investigations in this direction (for a 
review, see e.g. 0). All known thick br ane solutions can be classified into two large groups: static solutions and 
solutions depending on time (cosmological thick branes). From the mathematical point of view, a search for thick 
brane solutions implies the derivation of regular solutions of multidimensional field equations with a given material 
source. As a source of matter, various scalar fields are often used. The asymptotic geometry of spaceime outside a 
brane (the bulk geometry) can be either flat, de Sitter or anti-de Sitter. 

One of the main advantages of use of scalar fields in modeling thick branes is that on the one hand the models 
obtained are comparatively simple, and on the other hand scalar fields are fundamental by their nature. However, this 
does not exclude the possibility of creating brane models in using other fundamental fields. For example, it could be 
spinor fields. As we know, there are no so many solutions with spinor fields in general relativity. In particular, some 
cosmological solutions with spinor fields were found in Refs. [3]- 6] . In this connection, a search for other solutions 
would be interesting in itself. As applied to branes, in Q we considered a static model of a thick brane supported 
by a nonlinear spinor field. In this model, the five-dimensional regular solutions both with and without cosmological 
term were found. 

The solutions obtained in Q are completely static, when all physically measurable quantities (the metric function 
and the spinor field) are to be time- independent. The next step is to examine a model when a spinor field is time- 
dependent. In this case, since spinor field components are bilinear combinations such as ipip, one can seek a solution 
for the spinor field in the form ip(r,t) = /(r)exp 1 "', where f(r) is a column matrix which components depend only 
on the extra space coordinate r. The presence of such time-dependent components of the spinor field leads to the 
appearance of the off-diagonal components of the energy-momentum tensor. To avoid this, we add to the model one 
more spinor field. In this case one can choose two spinor ansatze in such a way that the off-diagonal components of 
the energy-momentum tensor of every spinor field would have the opposite signs and annihilate each other. 

The paper is organized as follows: in section [TTJ we derive the field equations for the model of a thick brane with 
two time-dependent spinor fields having an arbitrary potential energy, and with a cosmological term taken into 
account. Choosing the special ansatze for these fields, the off-diagonal components of the energy- momentum tensor 
are excluded. In section Hill an example of numerical calculations for a special form of the potential energy is presented. 
In section [TVl the trapping of a test scalar field to the brane is considered. 
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II. FIELD EQUATIONS 

We consider the five-dimensional gravitation with a nonlinear spinor field as a source of matter. The Lagrangian 
of the spinor field is 



ip flip - ipytp) - mipip + V(ip, tp). 



(1) 



We choose the potential V(tp,ip) satisfying the condition 

The corresponding five-dimensional Einstein and Dirac equations are 
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where a = 0,1,2,3,5 is the Lorentz index; A = 0,1,2,3,5 is the world index; e a A is the 5-bein; T a are the five- 
dimensional Dirac matrices in a flat Minkowski space; Dai[> = [pA — i w A ^afc) 4> is the covariant derivative of the 
spinor tp; r a b = i (r a Fb — rf,r a ); JJ^ = e A, y a DAtp; to, A are some parameters; A is the cosmological constant; r] a b is 
the five-dimensional covariant Minkowski metric. The energy-momentum tensor for the spinor field is taken according 
to the textbook H 
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where T A = e A T a are the five-dimensional Dirac matrices in a curved spacetime; g AB = e a A e B r] ab is the five- 
dimensional contravariant metric tensor; {p = ^T is the Dirac conjugated spinor; Da4> = 4> ( d a + TW A ab r a {, J with 



\p d a — Qa^P- Our definition of the energy-momentum tensor Q has the opposite sign comparing with [8[ in order 
to be consistent with the definitions for R a i> from Q. 

The five-dimensional Dirac matrices in a flat Minkowski space arc 
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where I 2X 2 is 2 x 2 unity matrix, and o\ are Pauli matrixes 
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We seek a wall-like solution of the system To do this let us choose the five-dimensional bulk metric in the 

form 



ds 2 = 2 (r) [e 2 ^dt 2 - dx 2 - dy 2 - dz 2 ) 



dr l 



(8) 



Next, we are going to consider a system consisting of two spinor fields ip\.2- The reason for this is to exclude 
off-diagonal components of the energy-momentum tensor. For the spinor fields, we use the following time-dependent 
ansatze 




and ip2 = e" 




(9) 



3 



Then the corresponding components of the energy-momentum tensor Q arc 



(T : 
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It is easily seen from expressions pip that the contributions from ipi and ip2 have the opposite signs for the off- 
diagonal components Tgg . Consequently, the combined energy- momentum tensor for two spinor fields has the diagonal 
components only. This allows the possibility of writing the Einstein-Dirac equations in the form 
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Using the redefinitions r/^fx — > r, oj^/h — > u, X/x — > A, K 3 a 4 — > a 4 , x 3 fo 4 — > b A , xA — > A, m^fn — > m and performing 
some manipulations with equations (| 14|) - (| 18|) . we obtain the following set of equations: 
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From equation pip , we see the following constraints on the boundary conditions 
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(19) 
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where the index zero refers to the quantities on the brane, i.e. when r = 0. 

The spinors ([9]) are chosen in the spinor representation, and in the standard representation the spinor components 
have to be either odd or even functions. Using this fact, henceforth we will assume that a = b a . 



III. NUMERICAL CALCULATIONS FOR A SPECIAL POTENTIAL V(tp,ip) 

In this section we consider a case when the potential V(ip,ip) from ([1]) is chosen in the form 
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FIG. 1: The profiles of the eigenfunctions a(r) (the solid line) FIG. 2: The profiles of <j>' {r) / <j>(r) and x{ r ) are shown by the 
and b(r) (the dashed line). solid and dashed lines, respectively. 



where [• • • ] denotes the antisymmetrization. For the potential given by (|25[) , we have 
V (V>i,2,V'i,2) = -4Aa 2 6 2 , 
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Using this potential, we can derive the following system of equations: 
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with the boundary conditions 



and the constraint 



a = b , (f>o = 1, (f>o = 0, xo = 0, Xo' = 
-4a 2 , (m + 2Aoq) 
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which follows from equation (|31|) . 

The numerical analysis shows that regular solutions of the set of equations (|2"9"l) - ([3"3l do exist only for an appropriate 
choice of the parameter uj. One can say that this system of equations describes a nonlinear eigenvalue problem for the 
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cigcnfunctions a(r), b(r). Then the statement of the problem is as follows: Given some values of ao, m and A which 
enter into equations (|2T))) - (j3"3")) and boundary conditions find the corresponding eigenvalue of u. We solved the 
system numerically using the NDSolve routine from Mathematica. As an example we used the following values of the 
parameters: ao = 0.9, m = 1.0, A = —0.25 which give us w ~ —1.45227015. The results of the numerical calculations 
are shown in Figs. [I] and [2] 



IV. TRAPPING OF MATTER 



The five-dimensional localized wall-like solutions obtained above can be used for a description of a brane only if it 
will be possible to show that various test matter fields can be confined on such wall. As an example of such field, let 
us consider here a test complex scalar filed rj with the Lagrangian 



L 



-d A r) d n- -m ry 77, 



where too is the mass of the test field. Using this Lagrangian, we find the following equation for the scalar field 

1 d 



(36) 



Here a/— 5 g is the determinant of the five-dimensional metric gAB, and n is a function of all coordinates r\ = r\(x A \ 
Taking into account that the canonically conjugate momenta = (E, ~p*) are integrals of motion, we seek a solution 
in the form 

t](x A ) = AT(r)exp(-ip Al x A1 ). 
Substituting this ansatz in (|36p . one can find the following equation for X(r) 



X' 
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we have 
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where the prime denotes differentiation with respect to r. To find a solution of this equation, let us determine the 
asymptotic behavior of the metric functions <f> and x- Taking into account the numerical results obtained in the 
previous section, one can see that asymptotically, as r — > 00, \ ~ * const. Then, using equation pip , one can easily 
find that </3 ~ exp(-y/— A/6r). Obviously, this expression is only valid for negative A. Correspondingly, equation (|37[) 
takes the following asymptotic form 
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with the asymptotically decaying solution 

Xoa ~ D exp 
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(38) 



where D is an integration constant. One can see from this solution that the test scalar field decreases exponentially 
fast that corresponds to the fact that this field is concentrated around the brane. 

As a necessary condition for the trapping of matter on the brane, one can require converging the field energy per 
unit 3-volume of the brane [loj . i.e., 
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and also the norm of the field rj should be finite 

oo 

HI 2 = I \^gv*Vdr= ( (f> i e x X 2 dr. 

J J —oo 

— oo 

Taking into account the asymptotic solution (|38p . one can see that both i^tot and H77H converge asymptotically. Then 
it becomes obvious from the above analysis that the localized solutions obtained in section Mil confine the test scalar 
field, and this indicates that such solutions can be interpreted as branc solutions. 

Summarizing the results, we have obtained the ^-symmetric thick brane stationary solutions supported by two 
nonlinear spinor fields in the presence of the five-dimensional cosmological A-term. We have shown that such solutions 
do exist for the special spinor ansatze given by (|9]). This ansatze allowed us to exclude the off-diagonal components 
of the energy-momentum tensor. It happens because the currents for the spinor fields are 

= W = ^^, (40) 



(^J 3 = -(^) 3 = ^— — , (41) 



a 2 -b 



and consequently 



J° = (J., i ) o + (J. 02 )°=2^^^j, (42) 
J" = (J^) 3 + (J^ 2 ) 3 = 0. (43) 

This means that the current J A along the z axis for the case of two spinor fields is equal to zero. This in turn leads 
to the vanishing of the off-diagonal components of the energy-momentum tensor. 

Using the ansatze ©, we performed the numerical calculations of the set of equations (f2"9l) - (|3"3"|) with a special 
choice of the potential V{i>,i^) in the form (|25[) . Equations (j2U)) - (j53")) constitute an eigenvalue problem for oj subject 
to boundary conditions (|34[) and constraint (1351) . The example of regular solutions of the above system is shown in 
Figs. Q]and[2] The asymptotic behavior of the obtained solutions corresponds to an anti-de Sitter spacetime (A < 0). 
Using the asymptotic solutions, in section Hvl it was shown that they can trap the test scalar field. It indicates that 
such solutions may be interpreted as brane solutions. 
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